from time import time
from scipy import *
from scipy.special import j0
import numpy
from sympy import *
import sympy

n = 11
ind1 = arange(n-1)
t0 = cos(pi*(ind1+1)/n)       
def f(l,j): return t0[l] - t0[j] 
print f(2,4)

x = sympy.Symbol('x')
f = sympy.diff(cos(x), x)

print type(f)
def f(x, y):
    return x + y

x = arange(10)
y = arange(10)

xx, yy = meshgrid(x, y)
print xx
print yy
print numpy.sin(xx + yy)

a = arange(1000)

print j0(a).sum()
print "enter parametrization:"

print eval("sin(3)")

t = time()
for i in range(1000):
    for j in range(1000):
        pass;

t2 = time()

print "time passed: ", t2-t